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Explicit expressions for all the transport coefficients have 
recently been found for a trapped Bose condensed gas at finite 
temperatures. These transport coefficients are used to de- 
fine the characteristic relaxation times, which determine the 
crossover between the mean-field coUisionless and the two- 
fiuid hydrodynamic regime. These relaxation times are eval- 
uated as a function of the position in the trap potential. We 
show that all the relaxation times are dominated by the colli- 
sions between the condensate and the non-condensate atoms, 
and are much smaller than the standard classical collision 
time used in most of the current literature. The 1998 MIT 
study of the collective modes at finite temperature is shown 
to have been well within the two-fluid hydrodynamic regime. 

PACS numbers: 03.75.Fi, 05.30.Jp, 67.40.Db 



Recently, the complete set of two-fluid hydrodynamic 
equations were derived starting from a microscopic model 
for a trapped Bose-condensed gas at finite temperatures 
[Q. These included all the transport processes which 
arise in a Bose superfluid, namely the thermal conduc- 
tivity, shear viscosity, and the four coefficients of sec- 
ond viscosity. Working within the Thmas-Fernii approx- 
imation, these two-fluid hydrodynamic equations can be 
written in the well-known Landau-Khalatnikov form plo] 
used to describe the hydrodynamic behavior of superfluid 
"^He. In the present paper, we use the explicit expressions 
for the transport coefficients obtained in Ref. M to de- 
fine three characteristic relaxation times associated with 
distinct relaxation processes due to collisions between 
atoms. As with a classical gas, relaxation times which 
are associated with transport coefficients play a crucial 
role in understanding the dynamics of a Bose-condensed 
gas. For example, they determine how fast thermal equi- 
librium is reached when the gas is perturbed, the rate of 
evaporative cooling in the Bose-condensed phase, and the 
characteristic frequency dividing the coUisionless from 
the collision-dominated hydrodynamic region. 

We recall that the transport coefficients describe the 
deviations of the non-condensate distribution function 
/(r,p,t) from the "local thermal equilibrium" form [Q. 
Thus it is natural that they describe the rate at which 
such deviations relax. Ref. ||l[ gives all the various trans- 
port coefficients explicitly as functions of the local fu- 
gacity, or, effectively, the local static condensate den- 
sity nco(r). However these expressions are quite complex 
and involve various integrals which must be done numer- 
ically. In the present paper, we evaluate the associated 
transport relaxation times as a function of position in 



the trap, for experimental parameters which have been 
used in recent experiments. In particular, we show that 
the pioneering MIT study B of collective modes at finite 
temperatures was well within the hydrodynamic region, 
contrary to what was originally estimated in that paper. 
In the absence of a quantitative theory of the collision 
times in a Bose-condensed gas, most of the BEG litera- 
ture up to now has used the Maxwell-Boltzmann elastic 
collision time |Q,^ 
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as an estimate of the relaxation rates in trapped Bose 
gases. In particular, ujTci = 1 has been used ^J 
to estimate the "crossover" between coUisionless modes 
{luTci 3> 1) and the collision-dominated hydrodynamic 
modes (wTd <^ 1). In (|l|), no(r) is the local den- 
sity of the thermally excited atoms, a = Sna^ is the 
Bose atomic cross-section in the s-wave approximation. 
V — y^8kBT/Trm is the mean value of the velocity of the 
thermal atoms in a classical gas (which is very close to 
that in an ideal Bose gas) . In the detailed solution of the 
Boltzmann equation for /(r, p, t) for a trapped Bose con- 
densed gas at finite temperatures |H,0,0|, an important 
distinction is made between collisions between thermal 
atoms (C22 collisions) and the collision integral involv- 
ing atoms being scattered into/out the condensate (C12 
collisions). In the classical limit, the C22 collisions give 
rise to the collision time given in (nl). C12 collisions only 
appear at temperatures below the Bose-Einstein transi- 
tion temperature Tbec a-nd are thus a unique feature of a 
Bose-condensed gas. Our detailed calculations in Ref. Q| 
and in this paper include both the C22 and C12 colli- 
sions. However, we find that the various transport relax- 
ation times for the thermal atoms are dominated by the 
C12 collisions in a trapped Bose gas because of the high 
condensate density. A key feature about trapped gases 
is that the condensate is strongly peaked (high density) 
at the center of the trap while the thermal cloud (non- 
condensate distribution) is always spread out over a much 
larger region. Even at T = 0.99TbeCj for example, the 
condensate density at the center of the trap is still com- 
parable to the non-condensate density. 

Our major finding in this paper is that the relaxation 
time T^ associated with the thermal conductivity is well 
approximated by a simple formula 
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as long as 
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(3) 



Moreover, numerical results show that all three charac- 
teristic relaxation times {t^^t^. and r^) arc of the same 
magnitude, within a factor of two. Thus we conclude 
that to a good first approximation, the relaxation time 
TBE(r) as defined in (g) can be used as a simple "univer- 
sal" collision time in trapped Bose gases, replacing the 
expression in ffl). One immediately sees that since (H) 
is valid in the center of the trap at effectively all tem- 
peratures below Tbec, we obtain rBE(r) •€. Tci(r). For 
modes concentrated in the condensate core region, it is 
TBE(r) which determines the crossover frequency between 
the collisionless and hydrodynamic frequency regions. In 
the thermal gas outside the condensate core, Td in @j 
gives a good estimate for the transport relaxation times. 
We now turn to the detailed calculations which lead 
to the above conclusions. Referring to Appendix B of 
Ref. ||l|, we first summarize the formal expressions for 
the various relaxation times in a trapped Bose gas: 



1. Thermal conductivity 
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where zq = e'^[''=°^'^°('')l is the local fugacity and gn{z) 
are the usual Bose-Einstein functions ^. The effective 
potential C/o(r) = C/cxt(i") + 2g[nco{v) -f no{r)] includes 
the external trap potential C/cxt(r) as well as the self- 
consistent Hartree-Fock mean field. The interaction be- 
tween atoms is treated by the usual low energy s-wave 
scattering length approximation, with g = AttH a/m. 
The thermal relaxation time Tk in (H) is given by 
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with 



Tk,12 = tbe- 
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where Td and tbe are defined in p) and (0) respec- 
tively. The equilibrium densities of the condensate and 
non-condensate are given by no(r) = (1/Aq)(73/2(2o) 
and nco(r) = [/Ltco - C^oxt(r)]/5 - 2no(r), where Aq = 
(2TTh^ /mkBTy^''^ is the thermal de Broglie wavelength. 



The dimensionless integrals I22 and 7^2 are defined in 
Ref. ||l[ and involve an integration over the C12 and 
C22 collision cross-sections (which involve products of 
the equilibrium Bose distributions describing the ther- 
mal atoms with an appropriate weighting function). 



2. Shear viscosity : 
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where the viscous relaxation time r^ is given by 
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Again, the dimensionless integrals /[j and /22 ^^^ given 
by somewhat complicated expressions involving products 
of the equilibrium Bose distribution functions |l| . 

3. Coefficients of second viscosity : 

From general considerations, Landau and Khalatnikov 
showed Q that in a Bose superfluid, there are four co- 
efficients of second viscosity. We find in our model that 
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where the relaxation time r^ describes how fast the chem- 
ical potentials of the condensate and non-condensate 
components come to equilibrium. 

The relaxation time r^, which arises entirely from the 
C12 collision integral in the kinetic equation, was first 
discussed in detail in Refs. R^. The hydrodynamic 
renormalization factor an in (E2) only involves the static 
thermodynamic functions of the gas, and is defined in 
Refs. [pP]. More explicitly, one finds 
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where T12 is a collision time associated with the C12 col- 
lisions nm- We note that in a Bose-condensed gas, the 
coefficients of second viscosity arise entirely from the C12 
collisions, in contrast to the k and 77 transport coeffi- 
cients, which also have a contribution from C22 collisions 
(which, however, is very small in a trapped gas). 
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FIG. 1. The Thomas-Fermi density profiles along the z 
axis of the condensate and the non-condensate at T = 1/iK 
(T ~ O.TTbec) with iV = 4 X 10'' atoms. 

From Eqs. (||) and (0), one sees that tbe/tk,i2 and 
Tci/tk.22 are given as universal functions of the local fu- 
gacity zo{r). It turns out that tbe/t,^,12 — 1 for zq ^ 0.3 
and Tci/rK,22 ~ 8/15 for < zq < I. Similarly, we see 
from Eqs. ( pX)| ) and ( pi] ) that TBEc/'''r;,i2 and rci/r^_22 are 
also universal functions of zq, but now have a more com- 
plicated dependence, increasing as zq increase. These 
analytical results are confirmed in Figs. 2 and 3. We em- 
phasize that while tbe in (0) gives a good estimate of 
Tk,,12 in (p), it is an empirical formula. 

We note that the transport relaxation rates involve in- 
tegrations which suppress the weight of the low energy 
collisions. As a result, there is no Bose enhancement 
from these processes and no divergence in 1/t„.,j when 
Zq — 6"^^""="'^'") approaches unity (at the edge of the con- 
densate or at Tbec)- This is in contrast with the collision 
times Ti2 and T22 defined as in Rcfs. [g||3], in which l/ri2 
and 1/t22 have an infrared divergence for zq — I due to 
the contribution from low momentum collisions. How- 
ever, note that 1/r^ in ( |l3|) has no divergence at zq — 1 

§■ 

As a concrete illustration of the behavior of these var- 
ious transport relaxation times, we consider the MIT 
study of collective modes at finite temperatures Q . The 
trap frequencies were lOr^jli: = LOy/2'iT ~ 230Hz and 
(jJz/^Ti — 18Hz, the scattering length a ~ 2.75nm for 
^^Na, and the total number of atoms was iV = 4 x 10^. 
Both the m — Q quadrupole oscillations and the out- 
of-phase dipole oscillations were studied in Ref. [^. In 
both cases, the condensate and the thermal cloud oscil- 
late along the 2-axis with the collective mode frequency 
of the order of the axial trap frequency lOz- It is therefore 
convenient to compare the various transport relaxation 
rates with the trap frequency ujz- 
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FIG. 2. The position dependence of the various dimen- 
sionless relaxation rates at T = l/j,K (T ~ O.TTbec )with 
A^ = 4 X 10 . The open circles show the approximate re- 
laxation rate I/wzTre using Eq. (Q) and the dashed line is 
the classical collision rate l/io^Tci using Eq. (|l|). 

In Fig. 1, we plot the spatial dependence of both the 
condensate and non-condensate densities along the axial 
z-direction at T ~ O.TTbec- In Fig. 2, we plot the spatial 
dependence of the various relaxation rates. We find that 
all the relaxation rates are much larger than the trap fre- 
quency LOz, and thus the collective oscillations with the 
frequency w ~ Wz are hydrodynamic modes. In Fig. 2, 
we also show the approximate relaxation rates using the 
simple formula given by (^ . This empirical expression is 
seen to be in excellent agreement with the thermal relax- 
ation time r„. For comparison, we also plot the classical 
elastic collision time given by (P. Clearly this classical 
collision rate is much smaller than all the transport relax- 
ation rates due to C12 collisions. Outside the central con- 
densate region, t^^ vanishes and t~^ and r~^ only have 
contributions from the C22 collisions. They can be well 
approximated by the classical gas results r„ ~ (15/8)rci 
and Tjj ~ (5/4)rci. To see the dependence on the total 
number of the atoms, we plot in Fig. 3 the relaxation 
times for N — 10^. In this case, all the relaxation rates 
are such that lu^t ^ 1 and thus the collective modes 
would be in the transition region between the collision- 
less and hydrodynamic domains. 

We might also note that the diffusive equilibrium re- 
laxation time Tp is always slightly smaller than the other 
relaxation times (see Figs. 2 and 3). As we have discussed 
in Ref. [p|M| , in principle one could have wr^ 3> 1 even in 
the hydrodynamic region where cuTf^^uTri <IC 1. In a uni- 
form Bose gas, this can occur in a fairy wide region near 
Tbec where r^ becomes very large pi . In this region, one 
expects to find a new zero-frequency relaxational mode 
|g|,||. However, the high condensate density in the cen- 
ter of a trap (relative to the low density of the thermal 
cloud) precludes one from observing this diffusive relax- 



ational mode - unless one works extremely close to Tbec- 
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FIG. 3. Plot of the dimensionless relaxation rates for 
N = l(f a.tT= 160nK (T ~ O.STbec)- 

One can evaluate the peak value of the approximate 
relaxation rate in y^ using the T = Thomas-Fermi 
expression for the peak density nco(O) at the center of 
the trap. This gives the useful formula: 
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where ano = {h/mOj)^''^, uj = (uj^^y^z)^ is the av- 
eraged frequency of the anisotropic trap, Nc is the to- 
tal number of the atoms in the condensate, and Tj^ = 
(hLu/kB)[N/({3)]^^^ is the transition temperature of an 
ideal trapped Bose gas. For a given value of N, this ex- 
pression can be used to estimate at what temperature we 
pass from the coUisionless to the hydrodynamic region for 
a given trap (wtbeCO) = !)• 

To be in the two-fluid hydrodynamic region (wtbe < 
1), one must work with a sample size with TV r^ 10^ 
atoms or more, when dealing with ^^Rb or ^^Na atoms. 
One can also make use of a highly anisotropic trap poten- 



tial UJx 



3> (jJz to have low-frequency collective mode 



Lu ^ LUz while keeping the high density due to the tight 
confinement in the radial direction. We note that if one is 
already close to the hydrodynamic domain above Tbec 
as estimated by using (|l|) ||§,|l^, one will be deeply in 
the hydrodynamic region (wtbe ^ 1) when the Bose- 
condensate forms. This is simply a result of the fact that 
in the condensate core region, tbe <C Tci- 

The scattering length of metastable ^He* |10| is esti- 
mated to be a ~ 16nm, much larger than in *^Rb or ^"^Na. 
Combined with the large sample size {N ~ 10^), the re- 
cent Bose condensates observed in *He* are well within 
the hydrodynamic domain even above Tbec ■ Finally, the 
adjustable value of a achievable near a Feshbach reso- 



nance (see Ref. [0) could be used to study collective 
mode frequencies as one smoothly goes from the mean- 
field (small a) to the collision-dominated hydrodynamic 
(large a) region. 

Detailed predictions for the collective mode frequen- 
cies in the two-fluid hydrodynamic region of superfluid 
trapped gases will be presented elsewhere, based on the 
theory given in Ref. [Q. In trapped gases, one has the 
interesting situation that while the thermal cloud in the 
condensate region may be well in the hydrodynamic re- 
gion, one eventually always enters the coUisionless region 
in the low density tail of the thermal cloud. This has not 
been included in the derivation of the two-fluid hydrody- 
namic equations given in Ref. ^. While this has little 
effect on the frequency of a condensate mode, it may be 
important in the damping of certain modes. Such effects 
have been included above Tbec by using a cutoff ||l2|,|l3[ . 
The main point, however, is that the damping in the 
hydrodynamic region will be directly proportional to an 
appropriate spatial average of the transport relaxation 
times introduced in this paper. 
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